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I. INTRODUCTION 

Transport in a two dimensional electron gas is affected 
by the spin-orbit coupling. The most common signature 
is the weak (anti)localization. It is a small correction 
to the conductance due to the interference of time re- 
versed trajectories. The sign of the correction depends 
on the presence or absence of the spin-orbit scattering. 
The correction is suppressed if a time reversal symmetry 
breaking magnetic field is presen1ii*2iii4. The spin-orbit 
term in the Hamiltonian has a form of a non-abelian vec- 
tor potential^. In the case of quantum dots, if the spin- 
orbit coupling strength is position independent, a gauge 
transformation can be done, which results in an effective 
Hamiltonian with reduced spin-orbit coupling, and a rich 
variety of symmetry classes^. If the spin-orbit coupling 
depends on the position the transformation can not be 
done any moreS^. As a consequence, the accessible sym- 
metry classes are the three standard classes of Dyson, 
classifying the systems according to the presence or ab- 
sence of time reversal and spin rotation symmetry. 

For quantum dots with chaotic dynamics random ma- 
trix theory gives a convenient way to describe the trans- 
port properties, provided that the electron transit time 
T erg is much shorter than the other time scales of the 
problem (mean dwell time Td w , spin-orbit time r so , mag- 
netic time tb, inverse level spacing) 9 . Constructing the 
appropriate RMT models for the crossover between the 
symmetry classes, the magnetic field and spin-orbit cou- 
pling dependence of the average conductance was calcu- 
lated in Refs. Ifjyla- The theoretical results are confirmed 
by numerical simulations— and they are in good agree- 
ment with the experiments 11-12 . 

If one would like to calculate the averages of other 
transport properties, such as the shot noise power, higher 
order cumulants of the distribution of the transmitted 
charge, or any other linear statistics, the density of trans- 
mission eigenvalues is needed 9 -. Jalabert et al>i£ gave the 
weak localization correction to the transmission eigen- 
value density for chaotic quantum dots belonging to 
Dyson's three symmetry classes. Our work extends this 
result to the crossover regime between these classes. We 



present a calculation of the dependence of the weak lo- 
calization correction to the transmission eigenvalue den- 
sity on spin-orbit coupling and perpendicular magnetic 
field. For the sake of simplicity we restrict our attention 
to symmetric dots, i.e. we assume that the two leads 
attached to the cavity support the same number of chan- 
nels N , and for technical reasons we consider the case of 
iV>l. 

In some other sense we complement the work done 
by Nazaro\ii^, who calculated the crossover behaviour 
of the weak localization correction to the transmission 
eigenvalue density for disordered samples in dimensions 
d = 1,2, 3. Here we give the results for d = correspond- 
ing to a quantum dot. 

The study of the transmission eigenvalue density is in- 
teresting not only because of the practical implications 
related to the linear statistics, but it is instructive by it- 
self as well, since it gives a deeper insight to the weak 
(anti)localization phenomenon. In the cases of Dyson's 
symmetry classes the weak localization correction to the 
transmission eigenvalue density is of the form of Dirac 
delta peaks at the endpoints of the spectrum^. Our ana- 
lytical, closed-form result shows that these peaks broaden 
in the crossover regime, but the correction still remains 
singular. Furthermore, similarly to higher dimensional 
cases of Ref. 0, it is possible to identify the peaks as 
originating from the singlet and triplet sectors of the 
Cooperon modes of conventional diagrammatic pertur- 
bation theory. Our result also enables us to study the 
transition from weak localization to weak antilocaliza- 
tion on the level of transmission eigenvalues. 

As applications, we calculate the weak localization cor- 
rection to the conductance, the shot noise power and 
the third cumulant of the distribution of the transmit- 
ted charge. In the case of the conductance we recover 
the result of Ref. 0, giving a verification of our calcula- 
tions. For the shot noise we find that for the symmetric 
cavities studied in this paper the weak localization cor- 
rection is absent in the full crossover regime. The third 
cumulant of the transmitted charge behaves the opposite 
way. It is "crossover induced" in the sense that the clas- 
sical contribution vanishesi^*i& and the weak localization 
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term is nonzero only in the crossover regime. 

The paper is organized as follows. In the next sec- 
tion we specify the systems under consideration and the 
model applied for the RMT description. We briefly sum- 
marize the formal definition and the practical importance 
of the density of transmission eigenvalues. In Sec. IIIII wc 
present our main result, the weak localization correction 
to the transmission eigenvalue density, and analyze its 
behaviour as a function of the degree of time reversal 
and spin rotation symmetry breaking. In Sec. UVI we ap- 
ply our result to the transport properties above. Finally 
we conclude in Sec. [V] 

II. DESCRIPTION OF THE SYSTEMS AND 
THE RMT MODEL 

Let us consider a chaotic quantum dot with two leads 
attached to it. We assume, that the number of propa- 
gating modes is the same for both leads. We choose the 
spin-orbit coupling to depend on the position to avoid 
the reduction of the coupling strength. The weak mag- 
netic field is perpendicular to the plane of the dot. The 
assumptions for the spin-orbit coupling and the magnetic 
field ensure that our system exhibits a crossover between 
Dyson's standard symmetry classes. 

The transmission eigenvalues are the eigenvalues of the 
matrix product t't'\ Denoting the number of modes in 
a lead with N, the transmission matrix t' , describing the 
transmission from lead 2 to lead 1, is an N X N matrix 
with quaternion elements. It is submatrix of S, the 27V x 
2N scattering matrix of the system: 

t' = W 1 SW 2 , (1) 

where W\ is an N x 2N matrix defined by {W\)ij = 1 
if i = j and otherwise, W 2 is a 2N x N matrix with 
{W 2 )ij = 1 if i = j + N and otherwise. The product 
t't'^ has 2N eigenvalues, where the factor two comes from 
the quaternion structure of the matrix elements. If the 
system is time reversal invariant, there are N twofold 
degenerate levels. 

We assume that the system can be described with ran- 
dom matrix theory, i.e. r crg <C Td w , tb, t so , where the 
magnetic time is related to the flux $ through the system 
as&ii 

i_ k /$y 

TB Terg \$0 ) 

where $0 is the flux quantum and n is a numerical factor 
of order unity. 

To give a statistical description of the crossover be- 
haviour of the transmission eigenvalues we need an RMT 
model for the scattering matrix in the crossover regime. 
This is provided by the "stub model"^ which was 
adapted for the system under consideration in Ref. 0. 
In this approach the S matrix is represented as 

S = PU{\- RU)- 1 P\ (2) 



with 

R = Q f rQ. 

In the above expression U is an M x M random unitary 
symmetric matrix taken from Dyson's circular orthogonal 
ensemble^ (COE) and r is a unitary matrix of size M — 
2N. The 2N x M matrix P and the (M - 2N) x M 
matrix Q are projection matrices with Py = <5j.j and 
Qij = $i+2N,j- The quaternion elements of the matrices 
U, P, and Q are all proportional to the 2x2 unit matrix 
12- The matrix r is given by 

r = c -TT H ', (3) 

where A is the mean level spacing of the dot. H' is an 
(M — 2N) dimensional quaternion matrix generating the 
perturbations to the dot Hamiltonian, 

H' = ixXt 2 + ia so {Ai<7 x + A 2 a y ). 

(4) 

Here Aj (j — 1,2) and X are real antisymmetric ma- 
trices of dimension M - 2N, with Tr (AiAj) = M 2 Sij 
and Tr (XX T ) = M 2 and <Ji are the Pauli matrices. 
The first term in Q describes the time reversal sym- 
metry breaking through the magnetic field. The second 
term, having a symplectic symmetry, corresponds to the 
Rashba and/or Dresselhaus terms in the case of position 
dependent spin-orbit couplingiS. The dimensionless pa- 
rameters x and related to the corresponding time 
scales as 

2 2ith 2 2ixti 
X = ^A' flso = ^A- 

At the end of the calculation the limit M — > 00 should 
be taken. 

The density of transmission eigenvalues is defined as 

P(T) =(£ i S(T-T i )) 

= -Ihm_ 0+ In^Tr^^i^)' , '" 

where the trace is taken over channel and spin indices. 
Having p(T) at hand we can calculate the ensemble av- 
erage of any linear statistics 9 

2N 

A = ^>(T,), 

i=i 

as 

(A) = J p(T)a{T)dT. 

Prominent examples for linear statistics are the con- 
ductance, the shot noise power, or the cumulants of the 
distribution of transmitted charge2i2fiiSl. The weak local- 
ization correction for the linear statistics can be obtained 
from the weak localization correction to the transmission 
eigenvalue density. 



III. RESULT AND DISCUSSION 

To find the density of transmission eigenvalues, one 
has to substitute the scattering matrix into the def- 
inition JSJl using lJT]l. expand the inverses and calculate 
the average with the help of the diagrammatic technique 
of Ref. |22 up to subleading order in the small parame- 
ter 1/N. The details of the calculation can be found in 
Appendix lAl The result is 

p(T) = Po (T) + 5p(T). 



With the 0(N) contribution 



2iV 



if < T < 1 



otherwise 



(6) 



we recover the known result of Refs. I9HTI2I24I The 

factor of two accounts for the spin, as a consequence of 
the known fact, that the O(N) term is insensitive to the 
spin-orbit interaction and the magnetic field. 

The main result of our paper is the closed-form, ana- 
lytical expression for the weak localization correction to 
the density of transmission eigenvalues, 



5p(T) 



tVt(I-T) 



E (- 

1=0,1 



rifa„) , r 2 (y n )-T 3 (y n ) 



) 



for < T < 1 and otherwise. The variables y„ 
ya = T and y\ = 1 — T. The expression for T m is 



(7) 



arc 



r m (y) 



7m (7m + 2 A) 
7^+4iV( 7m + A)y' 



with 



7i 



72 



2ai 



73 



Note that 5p(T) is antisymmetric with respect to the 
point T = 1/2. 

As a verification of J7J, we consider the limits of 
Dyson's symmetry classes. Labeling them in the usual 
way with the index f3, they correspond to x 2 , a 2 Q -C N for 
[3 = 1 i.e. systems with time reversal and spin rotation 
symmetry, x 2 ^S> N for (3 — 2, that is, systems where the 
time reversal symmetry is broken by the magnetic field 
and i 2 CjV<C a 2 a for f3 = 4, for time reversal invariant 
systems without spin rotation symmetry. In these limits 
expression Q recovers the known result 



Sp(T) 



2_-l 

' 4/3 



[S(T-0+)-S(T -1 + 0+)] (8) 



of Refs. Ii"ll.3l The factor of two in JSJ is due to the 
twofold degeneracy of transmission eigenvalues for [3 — 
1,4, mentioned in Sec ITU 

To get some insight how the weak (anti)localization 
peaks in (JSJ) emerge, let us have a closer look at 6p(T) as 




T 



FIG. 1: The weak localization weak antilocalization transi- 
tion for two values of the magnetic field, x/VN = 0.3 (a) and 
x/V~N = 1.3(b). We removed the parameter independent sin- 
gularities at the endpoints, 5p (T) = \/T(l - T)8p(T). The 
spin-orbit coupling strength is characterized by a ao /\ // ~N = 
0.0 , 0.3 , 0.4 , 10 for solid, dashed, dotted and dot-dashed 
lines respectively. 



it approaches the different limits. For j m <§; N the ex- 
pression for T m (y) is well approximated by a Lorentzian 
in the variable ^/y, 



r m (y) 



( lm /2N) 2 + (^y) 2 



In the opposite limit, j m 3> N, the functions T m become 
independent of y. Thus going to f3 = 1, T 2 and T 3 cancel, 
and Ti, together with the inverse square root prefactor 
evolves to the peaks in JSJ at the edge of the spectrum. 
Close to f3 = 2 the correction vanishes. Notice, that 
T3 is independent of the spin-orbit coupling parameters. 
Specially in zero magnetic field it always gives Dirac delta 
contributions at T — and at T = 1. Approaching (3 = 4, 
the contributions from T\, T2 disappear, thus the zero 
magnetic field peaks associated with T3 show up as the 
weak antilocalization correction. 

In the crossover regime the weak localization peaks 
in JHJ broaden, but the correction remains singular at 
the endpoints. If all the 7 m -s are finite, the singularity 
is present through a inverse square root factor propor- 
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tional to the O(N) contribution. The second, nonsingu- 
lar factor in Q determines the form of the weak local- 
ization correction as the function of the magnetic field 
and spin-orbit coupling. In the absence of the magnetic 
field (73 = 0) this picture is modified by the inclusion 
of the remanent Dirac delta contribution in place of the 
T3 term. On Fig. ^ we illustrated the transition from 
weak localization to weak antilocalization for two values 
of the magnetic field. The relevant regions are close to 
the endpoints, where the peaks of the pure symmetry 
cases (JSJ are located. Due to the antisymmetry of Sp(T) 
we plotted only the part around T = 0. 

The presence of spin-orbit coupling dependent and 
spin-orbit coupling independent contributions is analo- 
gous to the case of higher dimensional systems studied by 
Nazaro\ii^, where they correspond to contributions com- 
ing from spin-triplet and spin-singlet Cooperon modes. 
The situation is very similar in our case. The basic build- 
ing block of the diagrammatic expansion for Sp(T) is the 
combination TCT, with T = I2 (£> <Ji and 

C" 1 = M 1 2 <g> 1 2 - Tr R <g> R, 



where the function T is 



r(A,y) 



1 - 4iV 2 A 2 



(1 - 2NX) 2 +8NXy' 



The appearance of the spin-orbit coupling independent 
term is because of the decoupling of the singlet and triplet 
sectors of TCT. Going back to the discussion of the weak 
localization - weak antilocalization transition with 
in mind, we find that the weak localization peak is due 
to the triplet terms with eigenvalue \\±\, and the weak 
antilocalization peak comes from the singlet contribution. 

An other property of the higher dimensional cases 
that persists also for quantum dots is the breakdown of 
the perturbation theory for small magnetic fields near 
T = 0,1. More precisely, for fields x 2 < 1 the appli- 
cability condition dp -C p of the perturbation theory is 
violated in an interval of order 0(x/N 2 ) from the end- 
points of the spectrum. The conclusion is the same as 
in Ref. Q, namely, at small fluxes, for obtaining the de- 
tailed behaviour of the density near T = 0, 1 one has to 
treat the problem in a nonperturbative way. 



where () denotes quaternion complex conjugation and the 
tensor product is defined with a backwards multiplica- 
tion: 

(<7i (g) a 2 ){ai> ® a r ) = (ai<j v ) ® (oyer.,). (9) 

The trace in the second term is understood as 

(Trii® ~R) a p tl g = Rij,apRji. n s, 

where latin letters are channel indices, Greek letters re- 
fer to spin space and summation over repeated indices is 
implied. The very same structure emerges in the work 
of Brouwer et. ali£, where the authors identify C as the 
equivalent of the Cooperon in the conventional diagram- 
matic perturbation theory. In the limit M — > 00 it be- 
comes: 

CT 1 = 2{N + x 2 + a 2 so ){t2®t 2 )-a 2 so (a x ®(T x + <7y®ay). 

(10) 

If according to the multiplication rule we define the 
action of a matrix on a vector as (Av) a p = A aPia pv pa , 
the spin-singlet and spin-triplet basis turns out to be the 
eigenbasis of the matrix TCT with eigenvalues 

Aqo 1 = 2(7V + 73 ), Ali = 2(N + 7l ), A^ 1 = 2(N + 72 ). 

As in Ref. Il4l only the triplet eigenvalues depend on the 
spin-orbit coupling strength. The correction Q can be 
expressed as 

6 P (T) = 



-. 7= rp\ E (-1)" E r(A lm ,2M)-r(A o,2M) , 

(11) 



IV. APPLICATIONS 

Having obtained the 8p(T), let us see some applica- 
tions. In the following we assume zero temperature. First 
we compute the conductance. We find 



G = c ^S 1 a dT{p {T) 
1 



Sp(T))T 
, 1 



2( 7l +JV) ^ 4( 72 +Af) 4( 73 +W) 



where the second term represents the weak localization 
correction. It is another verification of Q, as we re- 
cover the corresponding result of Brouwer et. al. in 
Ref. [J Note that the correction is of the form SG oc 
Aoo-E™=_iAi m (Ref.S@). 

As a second application we consider the shot noise 
power. We get 



P 



2e 3 V 
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1 2e 3 V N 

dT( Pa (T) + Sp(T))T(l — T) = — - — — , 



that is, the 0(1) contribution from 5p(T) is absent. In 
the case of pure symmetry classes it is a known result, 
that the weak localization correction to the shot noise 
power vanishes if the number of modes is the same in 
both leads&ii It was shown in Ref. that this persists 
to the case of a p = 1 — * 2 transition too. Our result 
allows us to extend this prediction to the more general 
crossover interpolating between all of Dyson's symmetry 
classes. The reason behind the absence of the O(l) con- 
tribution is that T(l — T) is symmetric with respect to 
the point T = 1/2 and it is integrated with the antisym- 
metric density function 5p(T). 

To see an example, where the weak localization correc- 
tion is absent in the limit of pure symmetry classes, but 
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not in the crossover regime, let us take the third cumu- 
lant of the distribution of the transmitted charge. It is 
the opposite of the shot noise in the sense, that for cav- 
ities with leads supporting the same number of channels 
the O(N) term vanishesi^i^, thus the leading order of 
this quantity is determined by Sp. The third cumulant is 
proportional to 

h = I dT{p (T)+Sp(T))T(l-T)(l-2T). 
Jo 

The weak localization correction trivially vanishes in the 
pure symmetry case because of the factor T(l — T) and 
the Dirac delta functions in (JSJ . In the crossover regime 
we find 

f Ar f 7i(7i + 2iV) 72(72 + 27V) 73(73 + 27V) \ 
/3 I, 16(7i + AO 3 32( 72 + 7V)3 32(73 + J ' 

that is, the (ensemble average of the) third cumulant is 
"crossover induced" for a symmetric cavity. 

V. CONCLUSIONS 

We investigated the crossover behaviour of the weak 
localization correction to the density of the transmis- 
sion eigenvalues between Dyson's three symmetry classes 
f3 = 1,2,4 for a case of a chaotic cavity with symmet- 
ric leads. Using the stub model approach for the RMT 
description, with the help of the diagrammatic method 
of Brouwer and Beenakker—, we carried out a sublead- 
ing order calculation in the small parameter 1/N. Our 
main finding is a closed-form, analytical expression for 
the correction. 

We studied the weak localization - weak antilocal- 
ization transition in detail. We found that the weak 
(anti)localization peaks © of the case of pure symmetry 
classes broaden in the crossover regime, but the correc- 
tion still remains singular at the endpoint of the spec- 
trum. With our result Q at hand, we gave a quan- 
titative description of the broadening and the crossover 
from localization to antilocalization as the function of the 
magnetic field and spin-orbit coupling. 

We compared our results to the known cases of higher 
dimensionalities, and found strong similarities. First, our 
result also splits into spin-singlet and spin-triplet parts, 
with only the triplet contribution depending on the spin- 
orbit coupling. In the limits of pure symmetry classes, 
the weak localization peak comes from the triplet con- 
tribution, while the antilocalization peak is due to the 
singlet part. Second, we also find that for small mag- 
netic fields, the perturbation theory fails to describe the 
details of the density near the endpoints of the transmis- 
sion eigenvalue spectrum. 

We applied our results to the conductance, the shot 
noise power and the third cumulant of the distribution 
of the transmitted charge. The conductance served as a 
test for our calculations, we recovered the result of Ref. .7 



obtained in the framework of the same model. For the 
shot noise power we found that the weak localization cor- 
rection is absent in the full crossover, due to the symme- 
try of the transmission eigenvalue density. For the third 
cumulant we found opposite behaviour. It is crossover 
induced: the O(N) term is absent, and the 0(1) contri- 
bution is nonzero only in the crossover regime. 

Further directions of research could be to apply our 
result to obtain the weak localization correction to the 
full statistics of the transmitted charge. Another possi- 
bility would be to extend our calculations to the case of 
cavities with asymmetric leads. In that case, differently 
from the present results, we expect a nontrivial magnetic 
field and spin-orbit coupling dependence also for the shot 
noise power. 
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APPENDIX A: DETAILS OF THE 
CALCULATION 

In this appendix we give the details of the derivation 
of our main result (JJJ . We adapt a procedure of Brouwer 
and Beenakker— that removes the nested geometric series 
in JSJ , appearing due to the inverse in the expression (J2J 
for the S matrix. The price for this is the introduction of 
more complicated matrix structures. 

Let us introduce the 2M x 2M matrices 

(S o \ / c 2 \ (u \ 

F(2)= (f(z) F ] )> R= (0 fit)' 

( A1 ) 

where for S we use the representation (j3J) with P being 
M x M matrix 

Pij = 1 if i = j < 2N and otherwise, 

and the M X M matrices C\ and C 2 are 

(Ci)jj = 1 if i = j < N and otherwise 
C 2 = P-Ct. 

The Green functions F\{z) and F 2 (z) are defined as 

F(z) = C 1 (z-SC a S*C 1 )- 1 , (A2a) 
F'(z) = C 2 (z - tfCiSCn)- 1 . (A2b) 

The density of transmission eigenvalues can be obtained 
from F(z) as 

p{T) = -it- 1 lim lm.Tr{F(T + ie)). (A3) 
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The matrix Green function F(z) can be expressed as 
F(z) = 

{2z)- x E± (C ± C [1 - U(R ± Cz- 1 / 2 )] _1 UCz- 1 / 2 ) 
= (2Z)- 1 E± [C ± A±(F± - X±)B±] , 

(A4) 

withX± = R±Cz~ 1 / 2 andF± = X±(l - UXi)" 1 . We 
defined A± and B-t such that A±X± = C, X±B± = 
Cz" 1 / 2 . 

To get the ensemble average of F, one has to calculate 
the COE average of F±. In the following F± refers to 
this unitary average. It is related to the self energy E± 
through the Dyson equation 

F± =X±(1 + S±F±), (A5) 

We can express (F) directly through E± as 

(F) = (2Z)- 1 (C ± C ((1 - S ± X±)- 1 S±) Cz- 1 ' 2 ) 

(A6) 

First we calculate F± to leading order in ■h. To this 
order we have to consider the planar diagrams only. De- 
noting the resulting series as F± , for the self energy we 
find 

si 0) =E^(M° 3 ) \ (A7) 

n=l 

where the coefficients W n are given as 2 ^ 

w n = -N i - 2n (-ir- i ( 2n ~^). 

n \ n — 1 J 

The operator V acts on a 2M x 2M matrix A as 

(A n A 12 \ ( tvA 12 \ 

\A 2l A 22 ) ' ^trAsx J ' 

With the help of the generating function 

oo 

h(z) = Wnz^ 1 = — (y/M 2 + 4z - Af ) 

n=l 

we can write equation (|A7|I as 

= 

(PX±(1 - S^Xi)- 1 ) ft f (t>X±(1 - E^Xi)- 1 ) 2 ) . 
The solution is 

= ±(Vi-N^T)(; J). (A8) 
From i|A8fl it follows that 

27V 

TrF (z) = Tr^(z) = 




FIG. 2: Diagrams contributing to the explicit part JS^f of 
the weak localization correction to the self energy. 

from which we get the well known result (jSJ) for the den- 
sity of transmission eigenvalues. 

In accounting for the weak localization correction let 
us write the self energy as 

E± = 54 0) +6E±. (A9) 

It follows from i|A6|l . that F splits up too, 

F = F(°) + 5F, 

with SF containing the weak localization correction to 
the Green functions (|A2I) in its off-diagonal blocks. Up 
to first order in <$E±, after a little algebra we get 

*F = \ Q - ^) 5S ± C - ( A1 °) 

The contributions to the self energy correction <5£± 
come from the 0( A ^ a+i ) terms in the large-M expansion 
of E±. These can be sorted as 

oo 

ffi ± = <5E^ e) + ^VF n <J(PF ± ) 2n - 1 . (All) 

71=1 

The first term consists of diagrams with the outermost 
[/-cycle being non-planar (see Fig. |2J), and a term due 
to the sub-leading order in the large-M expansion of the 
cumulant coefficients 

oo 2 ^ 

<554 e) = ffiW + 8E&> + SW n (VF^) H ~ , (A12) 

n=l 

v ' 

with SW n = — ^~M2n — ■ Evaluating the diagrams of Fig. [21 
for 5£±' ) we find 

5E« = 

> (R^R + 2C 2 ) a , a " 
V [R^R + 2Ci) ct , ct aE aa ^ a i a <R^, a 

(A13) 



where ()* denotes the complex conjugate, Greek indices 
refer to spin space and we assumed summation for re- 
peated indices. Furthermore 



2v^T 



±- 



1 



Doing the summation in the third term in (|A12|) we get 



ST, 



(SW) 



b z - 1 / 1 
M z I 1 



(A16) 



-2Nb 2 TUT, G = a 2 T C^U T, (A14) 



with 



n = a 4 c- 



2\2\ 



{2Nb z ) 



The matrices T and C are defined as in sec. [rrn 

(2) 

The second term <5X± is 
6V (2) = ((* + S2 )Q n + 2 Sl Q 12 ) aa a>a ( J J 
where 

n _ , T ( a 2 G + b 2 E a 2 E + b 2 G 
y- M [ a 2 E + b 2 G a 2 G + b 2 E 1 N 



a 2 TrR®R* b 2 (M-2N) 
b 2 (M-2N) a 2 TrR®R*, 



where the trace is defined as in Sec. If III and 



si 



b ( z- 1 



3/2 



S2 



(A15) 



& 2 \ 
1 - 4z— si. 



The second term in (|AI 1 1) contains sub-leading order 
diagrams, that have planar outermost [/-cycles. Up to 
first order in 5T± 



S(VF ± ) 



2n-l 



V (FL 0) (5E ± F 



(«) 



VF 



(0) 



Putting everything together we see, that (|AI1|) is a (lin- 
ear) self-consistency equation for <5X±, which can be 
solved straightforwardly, if from i|A10|) we notice, that 
for the transmission eigenvalue density it is enough to 
get Tr 2 <5£±, where we denoted the spin-trace as Tr 2 . 
Substituting the solution in IjAlOfl in the lower left block 
we get the weak localization correction to F(z), from 
which using (|A3|) we arrive to the result J7J). 
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